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Abstract. We present the first in a series of papers that attempt to investigate the relation between binarity, magnetic activity, and 
chemical surface abundances of cool stars. In the current paper, we lay out and test two abundance analysis methods and apply 
them to two well-known, active, single stars, HD 27536 (G8IV-III) and HD 216803 (K5V), presenting photospheric fundamental 
parameters and abundances of Li, Al, Ca, Si, Sc, Ti, V, Cr, Fe, Co and Ni. The abundances from the two methods agree within 
the errors for all elements except calcium in HD 216803, which means that either method yields the same fundamental model 
parameters and the same abundances. Activity is described by the radiative loss in the Can H&K lines with respect to the 
bolometric luminosity, through the activity index i?HK- Binarity is established by very precise radial velocity (RV) measurements 
using HARPS spectra. The spectral line bisectors are examined for correlations between RV and bisector shape to distinguish 
between the effects of stellar activity and unseen companions. We show that HD 27536 exhibit RV variations mimicking the 
effect of a low-mass (m ~ 4Mj) companion in a relatively close (a ~ I AU) orbit. The variation is strongly correlated with 
the activity, and consistent with the known photometric period P = 306.9 d, demonstrating a remarkable coherence between 
RuK and the bisector shape, i.e. between the photosphere and the chromosphere. We discuss the complications involved in 
distinguishing between companion and activity induced RV variations. 

Key words, stars: abundances - binaries: general - stars: fundamental parameters - Stars: individual: HD27536 - Stars: 
individual: HD 216803 



1. Introduction 

There is currently still no single theory that adequately de- 
scribes magnetic activity on cool main-sequence stars with 
outer convective envelopes. Neither is there a solid description 
of its impact on stellar evolution. While the current scheme of 
dynamo models successfully account for the existence of mag- 
netic fields for stars at the age of the Sun, more quantitative 
explanations are lacking (see Strassmeier 2003. for a recent re- 
view). 

It has been established empirically ("Noves et alJ 1984 1) that 
the activity level, and hence the dynamo efficiency, is a mono- 
tonic function of the rotation period for stars with periods larger 
then a few days. It scales approximately with the Rossby num- 
ber - Pt^^, where Tc is the convective turnover time and 
P the rotational period. Hence, the activity is indirectly a func- 
tion of the spectral type. However, the boundary layer dynamo 
model is unable to account for the activity in the overactive M 
dwarfs, and also the reason for and the exact scaling with 'Rq is 
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* Based on observations collected at the La Silla Observatory, ESO 
(Chile), with the HARPS spectrograph at the ESO 3.6m telescope. 



not fu lly understood (e.g. ISchriiveilll996t ISchriiver & ZwaanI 
l200(i) . 

Although rotation is thought to be the necessary condition 
for magnetic activity, it is not clear to what extent binarity in- 
fluences the generation and morphology of magnetic fields and 
the coiTesponding chromospheric and coronal emission (e.g. 
Ibopp & Fekel 1977; Edwards 1983; Shkolniket al. 2003). The 
differential gravitational pull from a companion may cause a 
longitude- and latitude-dependent relationship between rota- 
tion rate and activity level, and may also contribute to an in- 
homogeneous chemical abundance by introducing tidal waves 
which affect the chemical stratification. No theoretical studies 
exist in the literature, but if such relationships exist, the models 
of the evolution of close binaries would then need to be recon- 
sidered. 

Another burning question is the link between chemical 
surface abundances and binarity. The recent observation that 
(solar-like) stars with planets are on average more metal rich 
than similar stars without planets (e.g. .Santos et al. 2001. 200^ 
l2005HFischer et alJ2003h is at least suggestive that this maybe 
also the case for binaries. Of particular interest is the fact 
that the Li abundance seem to differ between single stars and 
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stars which host giant planets jlsraeUan et aljEoC^ . and even 
between otherwise ide ntical twin-star binaries (|Marti'n et al] 
l2002t ball et alJ l2005'). possibly related to pre-MS accretion 
of planetesimals or engulfment of planets. 

Furthermor e, as instrumentation improves (i.e. HARPS - 
see iMavor et" al. 2003), more and more single stars turn out 
to have variable radial velocity (RV), due to the presence of 
low-mass companions, either low-mass stars, brown dwarfs or 
giant planets. It is thus becoming necessary and indeed possi- 
ble to deal properly with the effects of binarity, and to attempt 
to separate the classical dynamo process of internal differential 
rotation from the external contribution from close companions. 
In this discussion it is necessary to extend the concept of bina- 
rity to include all possible companions, including brown dwarfs 
and planets. 

We are currently conducting a search for true single stars 
among the known active stars, in order to study the activity- 
rotation relation in a sample that is not "polluted" by any type 
of binaries. Furthermore, we may hope to isolate the effects 
that binarity may have on stellar activity. Our sample consist 
of about 30 known active G-M stars, with known photometric 
variations attributed to rotational modulation of star spots, in- 
cluding both binaries and — supposedly — single stars. Most 
of our sample stars belong to the BY Draconis variables, which 
are defined from their variability due to changing star spot 
coverage (Bopp & FekeL 1977.). Several aspects need to be in- 
vestigated to properly analyze the sample, including (1) time- 
resolved RV monitoring and (2) line bisector analysis, (3) mon- 
itoring of activity indexes, (4) determination of fundamental 
atmospheric parameters, and (5) accurate photospheric abun- 
dance analysis. RV variations due to stellar acoustic oscilla- 
tions, which for our purposes is an additional noise source, are 
effectively filtered out by the long exposure times typically in- 
volved. 

In this paper we will investigate various tools available to 
perform reliable, accurate and fast abundance analysis. Since 
this will have to be done in detail for all stars in the sample, we 
are developing ways to perform this task in as automated a way 
as possible, without employing fully automated procedur es yet 
(e.g. minimum distance methods, see lAUende Prieto ^2004). al- 
though this is eventually the goal. Such automated tools will 
be of great utility in th e future, not least in view of upcoming 
facilities like STELLA JStrassmeier et alJ2004 . 

In this paper we present an analysis of the two well- 
studied active stars HD 27536 and HD 216803, deriving 
their fundamental atmospheric parameters and their elemen- 
tal abundances, as well as their activity indexes. So far only 
very few abundance analyses have been carried out for ac- 
tive stars, mostly due to the compl ications involved (see e.g. 
iKatz et alJl2003t iMorel et al.ll2003l for recent applications to 
RS CVn stars). We furthermore present a preliminary analy- 
sis of the time-resolved RV, bisector and activity index data for 
HD 27536. 

A future paper will deal with the full sample of stars ad- 
dressing the RV, activity, and line bisector variations and the 
statistical properties of the sample, presenting a discussion of 
the effects of binarity on the stellar activity. 
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Fig.l. The Can H&K hnes region of HD 27536 (top) and 
HD 216803 (bottom). HD 27536 is a G8IV-III giant while 
HD 216803 isaK5Vdwarf. 

Table 1. Log of the observations of HD 27536 and HD 216803. 
Exposure time is 400 s in all cases 



Star 


Date 


S/N @ 650nm 


HD 27536 


2004-10-01 


270 




2004-10-02 


420 




2004-10-31 


400 




2004-11-23 


245 




2005-02-12 


430 




2005-03-17 


230 


HD 216803 


2004-11-23 


230 



2. Observations 

2.1. The target stars 

HD 27536 (EKEri, HR1362 ) presents a special case. As shown 
by IStrassmeier et alJ lll999h it is a slowly rotating G8 giant 
with an over-active chromosphere possibly due to a fossil Ap 
type magnetic field, viewed at an inclination close to 90°. The 
star has a measured rotation period from long-term photome- 
try of 306.9+0.4 days. Its Li abundance is close to solar Yet 
its magnetic activity level is very high and comparable to some 
RS CVn binaries. 

HD 2168 03 (TW PsA , HR8721, G1879) is classified as a 
BY Dra star ( IBopd &Fekel 1977) . It is a K5 dwarf with a ro- 
tation period of 10.3 d ays but with a sm all projected rotational 
velocity of 2.6 kms"' dVoet et alJl983h . indicating an inclina- 
tio n of approxirnately 30°. The star has been studied recently 
bv lSantos et alJ(l200lh . who found atmospheric parameters and 
metallicity and placed it in a sample of comparison stars with- 
out planets. 

Both stars are thought to be single. 

2.2. New HARPS spectra 

Tabled show the log of the observations on HD 27536 and 
HD 216803. AU spectra were taken with the HARPS spec- 



T. H. Dall et al.: Binarity, activity and metallicity among late-type stars 



3 



trograph at the ESO 3.6m telescope at La Silla Observatory. 
Exposure times were 400 s in all cases, which ensures that 
stellar acoustic noise is effectively averaged out. In addition 
a high S/N solar spectrum was taken on the daylight sky. All 
spectra were acquired in the simultaneous ThAr mode, where 
a ThAr-lamp spectrum is recorded on the CCD alongside the 
object spectrum. The reduction process of HARPS spectra is 
fully automatic and extremely a ccurate due to the h igh intrin- 
sic stability of the spectrograph (Rup orecht et al . 2004). The fi- 
nal data products is a wavelength calibrated ID spectrum, with 
no normalization or flux calibration performed, and a cross- 
correlation function (CCF), computed for and averaged over 
all 72 spectral orders. The formal error on the wavelength cali- 
bration is ~ 0.02 ms"', while the total RV uncertainty is ~ 1 - 
2 m s ' for all spectra. The spectral resolution is R ~ 100, 000. 
HARPS spectra cover the range 3780 - 6910 A, except for a 
gap between 5304 - 5337 A where order 1 15 falls between the 
two CCDs. 

Both stars are clearly active, as evident from the Ca ii H & 
K lines, shown in Fig.^ 

3. Methodology 

The two methods for abundance analysis described in the 
following differ in their requirements onto the input spec- 
trum. While method 1 (Sect. lTTT l uses an automated procedure 
to fit the continuum and the lines simultaneously, method 2 
(Sect. l?^ uses an interactive procedure to normalize the spec- 
trum, followed by synthetic spectrum calculations. 

3.1. Abundance analysis: Method 1 

The abundance analysis consists of the following steps: mea- 
surement of the equivalent widths (EWs) over the full spectral 
region, determination of a first guess at the basic parameters of 
the star, calculation of the appropriate model atmosphere and, 
finally, the abundance analysis. The method is the same as em- 
ployed bv Dall et al.. (.2005..). adopted from the p rocedures of 
iMorel et alJll2003ll and lBmntt et alJ ll2002l kOOA 

The first step is the measure ment of the EWs, which is 
accomplished using DAOSPEC' dStetson & Pancino 2005. in 
preparation), which uses an iterative Gaussian fitting and sub- 
traction procedure to fit the lines and the effective continuum. 
The lines are identified using a list of lines from the VALD 
database jKupka et alJll999HPiskunov et alJll995h . where all 
lines deeper than 1 % of the continuum are included. Different 
line lists for different spectral types can be retrieved directly 
from the database. In Sect. 14.31 we discuss the choice of line 
list in more detail. For HD 27536 we used lines in the region 
5000 - 6800 A, while for HD 216803 the region 5500 - 6800 A 
was used. For bluer wavelengths the continuum determina- 
tion becomes uncertain, while the S/N decreases rapidly be- 
yond 6800 A because of the decreasing instrument efficiency. 
Next, an initial estimate of T^ff is found using the line depth 

' DAOSPEC has been written by P. B. Stetson for the Dominion 
Astrophysical Observatory of the Herzberg Institute of Astrophysics, 
National Research Council, Canada. 



ratios calibrated by iKovtvukh et all ll2003h . We approximate 
the line depths by using the EWs directly, assuming the same 
width for all lines. This is computationally much easier, since 
DAOSPEC does not provide us with the line depths. A more 
accurate determination of Tgg will be derived later in the pro- 
cess. When used on the spectrum of the Sun this procedure 
yields T^ff = 5793 + 37 K, hence using the EWs instead of line 
depths prove accurate enough for a first estimate. We next adopt 
log g and microturbulence parameter ^, for a canonical ZAMS 
star, and assume solar metallicity as our starting point. With 
this we then calculate t he initial model, using the ATLAS 9 
code adapted for Linux dKuruczll 19931 ISbordone et alJl2004 . 
With the measured EWs and the m odel, the abu ndances are 
calculated using the WIDTH9 code ( Ku ruczll99l modified for 
PC by V. Tsymbal), and compared line-by-line to Solar abun- 
dances, calculated from a high S/N solar spectrum. This last 
step is crucial to avoid problems due to uncertain ^/-values. 
The parameters of the model used to calculate the solar abun- 
dances are Tefr = 5778 K, \ogg - AAA, and ^, - 1.2 km s"'. 

Now the model parameters {T^s, \ogg, ^, [Fe/H], and 
[a/Fe]) are iteratively modified until consistency is reached, 
defined by the following criteria: (1) That there are no trends 
of Fe I abundance with EW, wavelength or excitation potential, 

(2) that the abundances derived from Fe i and Fe ii are the same, 

(3) that the derived metallicity and a-element abundances are 
consistent with the input model. 

The computational benefits of running these tools under 
Linux are immense: On a Pentium iii laptop the DAOSPEC 
computation of 800 EWs from the HARPS solar spectrum takes 
about three minutes, while the calculation of a new ATLAS9 
model takes two minutes. The abundance calculation using the 
modified WIDTH9 on 400 iron lines takes less than a minute. 



3.2. Abundance analysis: Method 2 

The second method is based on the calculation of synthetic 
spectra rather than just EW measurements. The software is 
called VWA which has a graphical user interface (GUI). The 
m ethod and the VWA softvy are have been described in detail 
by Bruntt et al. (2002, 2004) and the method compared to the 
classical EW method bv Bikmaev et al.. C200Z) . The abundance 
analysis relies on a t omic parameters from the VALD database 
jKupka et al."l99^: 'Piskuno v et alJ ll995^ and uses modified 
ATLAS9 atmospher ic mo d els fro m interpolation in the grid 
published by Heit er et al.l (l2002h . The least blended atomic 
lines are selected automatically by VWA but additional lines 
can be chosen manually. The abundance of each line is fitted 
iteratively by requiring that the EW of the observed and com- 
puted spectrum agree. The wavelength range used for match- 
ing of the EW is normally twice the FWHM of the line, but if 
a line is partially blended, non-symmetrical, or affected by a 
bad column in one of the wings, that part of the line may be 
excluded when calculating the EW. When all lines have been 
fitted the observed and synthetic spectra are compared. Offsets 
due to wrong placement of the continuum are easily identified 
by visual inspection in the GUI. 
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The VWA method requires that the spectra have been nor- 
malized. This is done by manually defining continuum win- 
dows which we fin d by comparing with a spectrum of the Sun 
llHinkle et al.l200(l) . and fitting a spline function through these 
points. 

3.3. Activity index 

We calculate the absolute emiss ion fluxes in the C a ii H & K 
lines, following the method by Linskv et al ] ill979ll : We inte- 
grate the emission flux, /k, in the K line between the Kiv and 
KiR points and we find /h in a similar way. This is then nor- 
malized by the flux /50 in the 50 A interval between 3925 and 
3975 A, and scaled to an absolute flux index: 

Tk = 5o|^r5o, (1) 

750 

where f^o is defined from Eq. 3 of lLinskv et alJ ( 1 19791) for V - 
R < 1.30; 

loglFso = 8.264 - 3.076(y - R). (2) 

Following lStrassmeier et all (I2OO0I) we use V - R colors com- 
puted via B - V colors from the color-color relation of Gray 
( I1992I) . since this minimize s the efifects of th e rotational mod- 
ulation of star spots ( Strassmeier et al.lll994l) . In order to cal- 
culate the activity index properly, we need to take into account 
the contribution of the photosphere. Hence the corrected flux 
indexes are (with a similar expression for the H line) 



[Fe-I/H] = 0.09 +/- O.OS 



(3) 



where is the index for a radiative equilibrium atmosphere 
without a chromosphere, as given in Linskv et al. ( 1979). We 
then find the 7?hk index as the total chromospheric radiative 
loss in the H and K lines in units of the bolometric luminosity: 



Rv 



o-T 



(4) 



eff 



where we use our derived Teff. Values of B -V are found via 
SIMBAD and checked for consistency with the derived T^s- 
All derived values are listed in Table|3 

4. Results 

4.1. Temperature, gravity, microturbulence 



Following the method outlined in Sect. 13.11 (Method 1) we 
start with the calculation of EWs directly from the pipeline 
produced ID spectrum, which results in 1095 measured and 
identified lines in HD 27536. Using the EW ratios instead of 
line depth ratios, we find Tetj = 5338 K ±60 K, and adopt this 
value plus \ogg = 4.40 as our initial guess for the model pa- 
rameters. After a few iterations we reach the model parameters 
Teff = 5240 K, logg = 3.55, ^, = 1.2 kms"', [Fe/H] = 0.09, 
[a/Fe] - 0.0, which fulfills the convergence criteria (see Fig.|3- 
For HD 216803 we aiTive at T^s = 4780 K, logg = 4.70, 
^, = 1.1 kms-', [Fe/H] = -0.05, [a/Fe] = 0.0. The uncer- 
tainties are estimated from the difference between the best fit 
model and the model corresponding to the Icr fits in Fig.|2l 



slope ■= 0;000131 



ope = -0.000991 



20 4D BO 80 

EW (mS) 



E.il.pol. (iV) 
[Fe-!!/H] = 0.10 +/- 



ope = 0.000002 



5000 
■do (S) 



Fig. 2. Method 1 : The iron diagnostic plots used in determining 
the atmospheric parameters of HD 27536, using the automati- 
cally selected Fe lines. Upper left plot and clockwise: Relative 
abundance of Fe i lines versus EW, same but versus excitation 
potential, relative abundance of Fe n Unes versus EW, Fe i abun- 
dance versus wavelength of the line. The dashed line is the 
mean value at zero slope while the dotted lines are the best 
linear fit and + Icr fits. 



Following Method 2 outlined in Sect. 13.21 we computed 
abundances for HD 27536, HD 216803 and the Sun. For the 
spectrum of the Sun, we find iron abundances of AA - H-0.09 
and -0.07 dex for Fe i and Fe 11 respectively, where AA is the 
abundance relative to standard value for the Sun, i.e. AA - 
log A^Fe /A^tot-(log A^Fe /Mot)©- The value for the solar abundance 
of Fe is (log A^Fe/Mot)© and is taken from iJrevesse & Sauvalj 
(1998) who found -4.54. Our measured abundances of Fei 
and Fe 11 are shown in Fig. |3] Open symbols are used for neu- 
tral lines and filled symbols are for ion lines. The abundances 
are plotted versus EW (left) and excitation potential (right). 
Note that only lines with excitation potential in the range 2 - 
5 eV and EW in the range 0-80 mA are plotted. In addition 
to the large difference between abundances of Fei and Fen 
we see a significant correlation of Fe i abundance with excita- 
tion potential, an effect also seen when using Method 1 . These 
discrepancies are explained by a wrong temperature structure 
in the ID-LTE ATLAS9 model for the Sun, as explained by 
Sh chukina & Trujillo Bueno (2 001). 

lAllende Prieto et all ll2004l) also found discrepant results 
for the Fe i/Fe 11 balance. The discrepancy they found increased 
when going to late type stars and they discussed whether spots 
in addition to the NLTE effects which are not included in the 
ID-LTE model atmospheres could explain this. In the present 
study we have used the ion-balance to modify the logg value 
of the model atmosphere, but this may not be the perfect ap- 
proach. However, as we show later, a strictly differential analy- 
sis, comparing the same lines between the stellar spectrum and 
the Solar spectrum, and adopting the abundance determined 
from the most abundant ionic species, will largely remove this 
problem. 
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Table 2. Parameters for the two stars. log/?HK and RV for HD 27536 are values from 2004-10-02, see Sect. l4.4l for details. See 
text for explanation of error estimates. 



Star 


T-eff [K] 


logs 


^, 


RV [kms-'] 


V sin I [km s '] 


[Fe/H] 


logifHK 




V -R 




HD 27536 


5240 + 35 


3.55 + 0.06 


1.2 + 0.1 


6.986 ± 0.002 


1.0 ±0.5 


-hO.09 


-4.33 ± 


0.12 


0.67 ± 


0.03 


HD 216803 


4780 ± 50 


4.70 ±0.10 


1.1 ±0.2 


7.228 ± 0.002 


1.5 ±0.5 


-0.05 


-4.62 ± 


0.21 


0.95 ± 


0.06 



Table 3. The abundances ([M/H]) derived by our two methods for HD 27536 and HD 216803 relative to the measured abundances 
in a Solar spectrum. The second column for each star lists [M/Fe], i.e. the abundance within the star of each element relative to 
iron. The errors given are the RMS around the mean. 









HD 27536 










HD 216803 








Method 1 




Method 2 




Method 1 




Method 2 






[M/H] 


[M/Fe] 


A' 


[M/H] 


[M/Fe] 


A' 


[M/H] 


[M/Fe] 


N 


[M/H] 


[M/Fe] 


N 


Al 


+0.18 


+0.09 


1 








-0.04 


+0.01 


1 








Ca 


+0.19(6) 


+0.10 


30 


+0.16(1) 


+0.07 


3 


-0.14(7) 


-0.09 


7 


+0.09(4) 


+0.08 


3 


Si 


+0.12(7) 


+0.03 


37 


+0.07(5) 


-0.02 


18 


-0.09(5) 


-0.04 


12 


-0.03(11) 


-0.04 


14 


Sci 




















-0.07 




1 


Sen 


+0.11(2) 


+0.02 


7 


+0.11 


+0.02 


2 


-0.01(9) 


+0.04 


5 


+0.03 


+0.02 


2 


Til 


+0.14(7) 


+0.05 


30 


+0.13(6) 


+0.04 


16 


+0.09(11) 


+0.14 


21 


+0.10(15) 


+0.09 


24 


Tin 








+0.07(11) 
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V 


+0.27(9) 


+0.18 


17 


+0.25(8) 


+0.16 


11 


+0.32(17) 


+0.37 


15 


+0.25(20) 


+0.24 


11 


Cr 


+0.13 


+0.04 


3 


+0.14(8) 


+0.05 


3 


-0.02(7) 


+0.03 


5 


+0.06(12) 


+0.05 


7 


Fei 


+0.09(6) 




377 


+0.09(6) 




73 


-0.05(7) 




144 


+0.01(8) 




78 


Pen 


+0.10(8) 




25 


+0.12(11) 




15 


-0.02(8) 




5 


+0.01 




2 


Co 


+0.08(9) 


-0.01 


9 


+0.00 


-0.09 


2 


+0.00(10) 


+0.05 


6 


+0.06(9) 


+0.05 


5 


Ni 


+0.07(8) 


-0.02 


122 


+0.06(5) 


-0.03 


21 


-0.08(7) 


-0.03 


49 


-0.03(9) 


-0.04 


23 



An alternative approach to the determination of Teif is the 
use of photometric color indices. Using the calibrations of 
lAlonso et al.l ('l999') for HD 27536 and the lAlonso et 111(119961) 
cal ibrations for H D 216803 with the Stromgren indices found 
by lOlsenI ll 19931) and B - V colors from SIMBAD, we find 
T^sib -y) ^ 5035 K, T^^iB - V) = 4998 K for HD 27536 
andTeifC^-y) = 4491 K, T^siB-V) = 4438 K for HD 216803, 
in both cases assuming Solar metallicity. The uncertainties 
are around 120 - 140 K for all values. Another tempera- 
ture indicator is the Stromgren ji index. However, no yS value 
is available for HD 27536. Instead, it was inferred from the 
spect ral type (/? = 2.54 + 0.01) using the Stromgren indices 
from lOlsenI (Il993l). We then use the TempLogG software^ 
(see e~ lKupka & Brunttll200lh which use the calibrations of 
iNapiwotzki et alJ lll993ir The result is logg = 3.1 + 0.7 and 
Teff - 5070 ± 130 K. The interstellar reddening was found to 
be E{b-y) = 0.045 + 0.008. The parameters of HD 216803 fall 
outside the valid range of the TempLogG calculations. The dis- 
crepancy between color-derived and ion-balance derived Teff 
points to problems with the model atmospheres, the most ob- 
vious on e being the lack of a chromosphere. As shown by 
IJVIorel et a l. (2003J, the use of ATLAS9 models without taking 
into account the chromosphere, results in an overestimation of 
Teff and logg. Their investigation was based on the merging 
of an empirical chromosphere model and an ATLAS9 model, 
and although crude their results showed that the added chro- 
mospheric heating required a photospheric temperature several 
hundred K cooler. They also concluded that the derived abun- 
dance ratios were not very sensitive to the addition of the chro- 



^ e.g. ,http://ams.astro.univie.ac.at/templogg/| 



mosphere. Lacking proper chromosphere models we have em- 
ployed the ATLAS9 models without modifications for consis- 
tency with past and future studies. 

To decrease the systematic errors from the NLTE effects 
and to reduce errors caused by wrong log gf values we com- 
puted abundances in HD 27536 and HD 216803 relative to the 
abundances found for the same lines in the Sun. This causes a 
decrease in the RMS of the abundance of Fe i by 40% (compare 
Figs.|3]and|3. We only use Unes with EW less than 90 mA and 
excitation potential in the range to 5 eV. The results of Fe for 
the best models of Methods 1 and 2 are shown in Figs.|2|and|4] 
respectively. 

4.2. Chemical abundances 

Abundances derived with our two methods are listed in Tabled 
and showed in Fig.|5] The errors quoted in the table are RMS 
around the mean abundance where a sufficient number of lines 
were used. Additional errors are introduced by uncertainties 
in the atmospheric parameters. Given the uncertainties on T^ff, 
logg and we estimate an additional error of 0.05 dex which 
is to be added quadratically to the RMS errors. Additional error 
sources are due to NLTE effects, which are known to aff'ect a 
few elements, like Fe (di scussed in Se ct. 14. 1> . V, Mn and to a 
lesser extent Ti ( Bodagheeetal ]l2003l) . For the remaining ele- 
ments the effects are unknown or uncertain. EW measurement 
errors are likely significant for elements with only a few mea- 
sured lines i.e., Al, Sc and Cr. Small systematic errors could 
be introduced by wrong continuum placement, but since we 
measure all abundances relative to the Sun, this contribution is 
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Fig.5. Relative abundances of HD 27536 and HD216803 from 
Method 1 (*) and Method 2 (o). Note the overall good agree- 
ment, the deviating vanadium abundance and the discrepancy 
for calcium in HD 216803. 



Fig. 3. Method 2: Absolute Fe abundance found for the solar 
spectrum (top), HD 27536 (middle) and HD 216803 (bottom). 
Open symbols are used for neutral iron lines and closed sym- 
bols are ionized iron lines. In the left panels abundances are 
plotted versus EW while in the right panels abundances are 
plotted versus excitation potential. The parameters of the ap- 
plied model is given in the insert boxes. 
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Fig. 4. Method 2: Relative (to Solar) Fe abundance for 
HD 27536 (top row) and HD 216803 (bottom row) plotted ver- 
sus EW (left column) and excitation potential (right column). 
Same models and symbols as Fig.|3] 

negligible. On the other hand, random errors from the contin- 
uum treatment will surely affect weak lines (e.g. Ti ii and some 
lines of Fen). 

The abundances of the main elements for several models 
are shown in Fig. |6l The left panel is for HD 27536 and the 
right panel is for HD 216803. The best model is labeled "F" 
while models with lower (higher) values of log g and tempera- 
ture are labeled "mG/mT" ("pG/pT"). The detailed parameters 



of each model is given below each column. For example the 
"pT" model of HD 216803 has AT = 100 K, Alogg = 0.0 and 
A^, =0.1 kms"' which is relative to the reference model (i.e. 
model "F") at Teif = 4800 K, logg = 4.8 and ^, = 1.5 kms-\ 
thus the temperature of model "pT" is 100 K higher and the 
microturbulence is 0. 1 km s"' higher than this. If there is a sig- 
nificant correlation (fitting a straight line) of abundance and 
excitation potential, a triangle symbol is used instead of a filled 
circle in Fig.|6l(e.g. Fe i for the "mT" model). Open circle sym- 
bols are used if less than five lines were used in the abundance 
determination (e.g. for Sc in both stars). 

A summary of the derived abundances of our preferred 
models "F" for HD 27536 and HD 216803 are given in Table|3l 

4.2.1. Notes on individual abundances 

In the following, we discuss the determination of abundances 
of individual elements and the error estimates. 

The lithium abundance For measuring the Li resonance 
line at /l/l 6707.8, the automatic deblending performed by 
DAOSPEC proved inadequate, and instead we used the IRAF^ 
task splot for the deblending as well as a spectral synthesis 
using VWA. The line is actually a dublet split by fine struc- 
ture. The major contribution comes from the more abundant ^Li 
isotope, with blends from weak Fe i, V i lines, which is taken 
care of with the VWA fitting, an example of which is showed 
in Fig. The equivalent widths measured are 15.5 mA and 
28.7 mA for HD 27536 and HD 216803, respectively. 

The Li lines in the solar spectrum could not be reliably mea- 
sured, so we are using directly the atomic parameters for the 
lines. The gf values of the Li lines are among the most accu- 

IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation, U.S.A. 
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Fig. 6. Relative abundances found for different models of HD 27536 (left panel) and HD 216803 by method 2. The preferred 
model is labeled F while models with lower (higher) \ogg are labeled mG (pG) and models with lower (higher) T^ff are labeled 
mT (pT). The differences with respect to model F are listed below each model. At the bottom of each panel the T^ff and \ogg of 
model F are given ("Reference"). Abundance values ai'e relative to the Sun. 
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Fig. 7. The lithium line region in HD 27536 with fits using three 
different abundances (-10.94, -10.84, -10.74). Blending with 
Fe I and V i, but not with CN, has been taken into account. 



The ff-elements The abundances of the electron donor ele- 
ments Ca, Si and Ti merit separate comments. In active stars 
these elements ar e often enhanced with respect to iron (e.g. 
iMorel et alJ2Q03l) . hence the enhanced opacity due the H ab- 
sorption need to be accounted for in the model calculations. 
While [Si/Fe] and [Ti/Fe] generally agree between Methods 

1 and 2, there is some discrepancy concerning calcium for 
HD 216803, possibly due to the abundance using Method 

2 being based on only three lines. Taking an unweighted 
mean of [Si/Fe], [Ti/Fe] and [Ca/Fe] for HD 27536 we 
get [a/Fe] = -h0.06 (-h0.03) from Method 1 (Method 2). 
Excluding calcium, we find similarly [a/Fe] - +0.05 (+0.03) 
for HD 2 16803. Within the uncertainties, this result is compat- 
ible with no (T-enhancement for both stars, and we did indeed 
use [a/Fe] = 0.0 in the model calculations. 



rate for any stellar lines iSmith et aLlliooi) . and can safely be 
used in an absolute manner. The result is A(Li) = -10.84 + 0.10 
and -1 1.19 + 0.13 for HD 27536 and HD 216803, respectively, 
the errors coming primarily from uncertainty in the continuum 
placement and the influence from the weak Fe i, V i, and CN 
blends. Identical results are obtained with splot and VWA. 

The contribution from CN to the total EW has not 
been taken into account, since we cannot establish reli- 
able carbon and nitrogen abundances. Only the two Ci 
/l/l 5052. 167, 5380. 337A lines can be measured, and they give 
deviating results. Assuming [C/Fe] = at most, we expect a 
contribution from CN of up to a few mA, resulting in a possi- 
ble overestimation of the hthium abundance by 0. 1-0.2 dex. 



Chromium Using Method 1 three Cr i and three Cr ii lines were 
used for HD 27536. The Cri lines yielded +0.13 with a spread 
of only 0.01 dex while the Crn lines gave +0.05, again with a 
spread of only 0.01 dex. Assuming this difference to be real and 
not just a statistical effect, we believe the difference to be due 
to shortcomings of the model atmosphere, rather than reflecting 
errors in log g. 

Vanadium For both stars we find very high abundances of 
vanadium ([V/Fe]) using either method. Rather than a true 
overabundance, this may actually reflect temperature inhomo- 
geneities on the stellar surface, indicating substantial spot cov- 
erage (e.g. iGrav & .lohanson .1991: .Strassmeier & SchordanI 
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Fig. 8. Same as Fig.[5| but using the manually selected line list. 



120001) ■ Vanad ium lines are know n to be strongly enhanced 
m sunspots ("Wa llace et al.l Il999^ since they have a very 
strong temperature dependence. Alternatively, the cause may 
be NLTE effects, but we do not find any trends of abundance 
versus EW, which would indicate this. 



4.3. Line list comparisons 

Several studies of stellar atmospheric abundances have em- 
ployed different line lists, all carefully crafted to suit the anal- 
ysis for which they are intended and this work is no excep- 
tion. Our line list is selected directly from the VALD database 
based on the quality of the spectrum under investigation, i.e. 
the list will depend on the S/N of the spectrum and the amount 
of blending, and hence be a function of spectral type and rota- 
tional velocity. We have tested two different selection schemes: 
(1) From the full line list we select manually the best non- 
blended lines from a comparison with the spectrum under in- 
vestigation, and (2) we use the full VALD line list and subse- 
quently reject outlying abundance values. 

The two schemes yield very similar results when used to 
find the stellar parameters for HD 27536. In scheme 2 we use 
346 Fei and 31 Fen lines after rejection iteration, while for 
scheme 1 we selected 205 Fei and 15 Feii lines. The differ- 
ences in atmospheric parameters for HD 27536 are AT^ff = 
10 K, Alog^ = 0.00, A^, - 0.1 kms ' between the two line 
lists. In Figs.|2and|Slwe show the diagnostic plots of the two 
different line lists, both using Method 1 for the analysis. The 
derived abundances differ by 0.02 dex at most, with most ele- 
ments differing less than 0.01 dex. Hence we can conclude that 
using all possible lines, relying on the statistical accuracy of a 
sufficiently large number of lines works as well as a manually 
selected set of lines. 

While the two methods prove to be internally consistent, 
we had some concern about how they compare with each other. 
In Fig.|5]we show the EWs measured using DAOSPEC (EWi) 
versus the EWs of the same fines using VWA (EW2), and as 



A A** 
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Fig. 9. The EWs as measured for 110 Fei lines using 
DAOSPEC (EWi) and VWA (EW2). The dashed fine is the 
EWi = EW2 relation. As can be seen, DAOSPEC EWs are sys- 
tematically lower than the synthesized. 
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Fig. 10. Iron diagnostic plot using the EWs found by Method 
2 with the model found by Method 1. See Fig.|2]for legend. 



can be seen, there is a systematic offset between the EWs of 
the two methods. 

In order to investigate the effects this offset may have on 
the model parameters, we calculated the abundances using the 
VWA EWs with Method 1. In Fig.^Jwe show the Fe diagnos- 
tic plot for the VWA EWs with the model found by Method 1 . 

As can be seen, the agreement is exceUent, and also the 
abundances are all the same within 0.02 dex, which means that 
either method of EW calculation leads to the same fundamental 
model parameters and the same abundances — provided they 
are properly calibrated on the spectrum of the Sun. We refer 
the interested reader to the DAOSPEC manual'* for an excellent 
discussion of continuum fitting and the problems involved. 

We necessarily also calculated the absolute iron abundance 
of the Sun, again using Method 2 EWs with Method 1; the 



* Can be obtained fro m_http://cadcwww.dao.nrc.ca/stetson/daospec/| 
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results of which is A(Fe)0 = -4.47 + 0.13, while using the 
DAOSPEC EWs we find A(Fe)0 = -4.53 + 0.14, using ex- 
actly the same set of Fe i U nes. The solar iron abundance from 
iGrevesse & Sauvall (Il998l) is A(Fe)0 - -4.54, hence it is not 
possible from this result to determine which EW determination 
is the most correct (if indeed such a concept is meaningful), 
since both methods give consistent results within the errors. 

4.4. Radial velocities, bisectors and activity indexes 

The activity indexes for HD 27536 and HD 216803 are Hsted 
in Table|5] The main error source is the errors onV - R, while 
the error from measuring the fluxes in the spect r a are around 
0.02 on log RuK for both stars. I Strassmeier et alJ lll990l) found 
logT^HK = -3.85 for HD 27536. Taking into account the dif- 
ferences i n derived TeM and V - R, the present result and the 
findings of lStrassmeier et aD ( ll99 0) are similar: using their pa- 
rameters, our values of log^HK would shift by -1-0.21. 

For HD 27536 we have acquired six RV measurements with 
irregular intervals (see Table These are shown in Fig. ^2 
along with an orbital fit using the known 306.9 d photomet- 
ric period. While the measurements are consistent with the pe- 
riod, it is clear that more measurements are needed. To inves- 
tigate whether this variation may be due to a low-mass com- 
panion, we rneasure d the line bisectors using the procedure of 
IOuelozetal.l ll200l'). We measured the bisector inverse slope 
(BIS) of the mean CCF of the spectra, and compared with the 
measured RV to check for correlations (Fig.^J upper panel). 
As can be seen there seem to be a correlation between the two 
quantities, indicating that the variation is internal to the star, 
probably due to activity i.e. spots moving across the stellar 
disk. Another indication of this comes from the activity index 
Rhk, which is shown in the middle panel of Fig. [2] as a func- 
tion of RV. This shows exactly the same behavior as the BIS, 
strenghtening the interpretation that the RV variation is due to 
activity variations. The correlation between the photospheric 
variation (BIS) and chromospheric (Rhk) activity seem to be 
very strict (Fig.[T21 lower panel), in accordance with expecta- 
tions (e.g. Messina et al. 2001). 

Obviously, part of the variation may be due to an unseen 
companion, hiding in the activity induced jitter Unfortunately, 
our measurements ar e too sparse to allow a definite conclusion. 
It is worth noting that lOueloz et alJ ( l200ll) found a negative cor- 
relation between BIS and RV in the case of HD 166435, which 
allowed the authors to conclude that the RV variation was in- 
duced by activity. On the other hand, we find a posit ive correla- 
tion between BIS and RV, as did Santos et al. '(2(1(1:1 in the case 
of the visual binary HD 41004, which was ascrib ed to a brown 
dwarf orbiting the fainter B component. Later, IZucker et alJ 
(12004') identified the RV signature of a giant planet around the 
A component. Hence, in the case of HD 27536 we cannot ex- 
clude the possibility of companion induced RV, activity and 
BIS variations, but more observations are evidently needed. 
Even so, if we assume that the activity is indeed induced by an 
orbiting companion interacting magnetically with the star, we 
would expect a RV "jitter" of 10-20 m s ' based on the calibra- 
tions of iSantos et al.. l2000i) . i.e. consistent with the observed 
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Fig. 11. RV variations measured in HD 27536. The variation 
is consistent with the known 306.9 day period, but obviously 
long-term monitoring is necessary. The dotted curve is a binary 
orbit fit with the known photometric period. 

variation, although their calibration strictly speaking only is 
valid for dwarfs. To disentangle any companion-induced pull 
from the activity jitter would thus require several years of mon- 
itoring, since the period is so long and since the orbital pull and 
the induced activity may have the same period. Assuming the 
period of the hypothetical companion to be equal to the photo- 
metric period of 306.9 d, and adopting an inclination of 90°, we 
find a mass of ~ 3.9 My orbiting in an eccentric (e ~ 0.5) orbit 
with semi-major axis ~ 1.2 AU. This is under the assumption 
that all the RV variation is due to an unseen companion. More 
likely, a companion would induce a smaller effect on top of the 
(larger) activity induced jitter. Note also, that the photometri c 
variation is sinusoidal, as shown by IStrassmeier et alJ ( Il999l) . 
which is hard to reconcile with the eccentric "orbit" suggested 
by the RV data. 

5. Conclusions 

In this paper we have presented a pilot study of active stars, 
presenting and discussing the tools to be used on a larger sam- 
ple. 

We have presented results on the two photometrically vari- 
able, late-type active stars HD 27536 and HD 216803, deriving 
their fundamental atmospheric parameters and photospheric el- 
emental abundances from high-resolution spectra. We accom- 
plished this task using two different methods, which produced 
results in good agreement with each other, giving credibility 
to the derived parameters and abundances. Furthermore, we 
have showed that, despite the obvious differences in derived 
EW between the two methods, and despite the differences in 
line lists employed, we derive essentially the same results. Both 
methods have a GUI, which greatly increases user-friendliness. 
Especially Method 1 has been developed with a view to au- 
tomation, and the ultimate goal is to be able to determine at- 
mospheric parameters and abundances for a great number of 
spectra without any human intervention. 

While the agreement between the two methods is very 
tight for HD 27536, larger differences resulted for HD 216803 
(A[M/H]~ 0.06). Even so, the abundances relative to iron, i.e. 
[M/Fe], show near perfect agreement between the two meth- 
ods, except for calcium, which may be due to poor statistics. 
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Fig. 12. Upper panel: Line bisector inverse slope (BIS; see text 

for definitions) vs RV. Middle panel: Activity index log/?HK 
versus RV. Lower panel: log/?HK versus BIS. There seem to be 
a clear correlation between activity level and RV, and between 
the BIS and the activity, pointing to changing activity as the 
cause of the RV variations. The error-bars on log^HK are only 
measurement errors i.e., not including uncertainties \nV - R 
and Teff , since these would affect all measurements equally. 



Our study demonstrates the importance of determining rel- 
ative abundances with respect to the Sun, which minimizes any 
errors introduced by uncertain gf values and uncertainties in 
continuum treatment. On the other hand, the presence of sur- 
face spots may influence the determination of logg through 
the ion balance, wherefore it may be better to use diff"erent ap- 
proaches, or select lines based on their sensitivity to spots. 

Due to the high intrinsic stability of the HARPS spectro- 
graph, a detailed analysis of the RV and bisector variations was 
possible, using the limited data set of HD 27536. We find a 
RV variation, which is consistent with the known photomet- 
ric period, although our coverage is not long enough to judge 
whether this is indeed the true period, or if the variations are 
periodic at aU. We find a remarkable correlation between the bi- 
sector inverse slope and the activity index log/?HK, which both 
vary in phase with the RV, reveaUng a very tight correlation 
between the photospheric and chromospheric activity. The fact 
that we find a positive correlation may hint at the existence of 



an unseen companion, but our phase coverage is not yet good 
enough to explore this possibility. 
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